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Abstract 

Describing and measuring the elastic properties of cellular materials such 
as honeycombs and foams can be a difficult problem when the cell structure 
is disordered. This paper suggests that tracking the flow of forces through 
the material can help in visualizing and understanding how the geometry 
of the cell structure affects the elastic behaviour. The mean strain tensor 
for a sample of material can be calculated by summing over the force paths, 
weighted by the strengths of the paths. This method emphasizes the paths 
with the greatest stress, which can have the most dynamic effect. The path 
averaging technique reproduces previous expressions for the Poisson's ratio 
of regular honeycombs, but easily extends to disordered honeycombs and 
foams. 
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1. Introduction 

Strain is a measure of how a material deforms when subject to stresses. 
For homogeneous isotropic materials the relation between stress and strain 
can be parameterized by elastic moduli such as the Young's modulus and the 
Poisson's ratio. When stretched in one direction most materials deform in 
the transverse directions. Poisson's ratio measures the relative amounts of 
strain in the parallel and transverse directions and is defined as positive for 
most materials where stretching in one direction causes contraction in the 
other two. 

Cellular materials such as honeycombs and foams can exhibit auxetic 
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(Evans et al. 1991) elastic properties, in particular negative Poisson's ra- 



tios (Lakes, 1987 Evans, 1989 Masters and Evans, 1996; Gibson and Ashby 



1997 Grima et al. 2006 Scarpa et al. 2008). For example, detailed model- 



ing of regular 2 dimensional molecular honeycombs (Evans et al., 1994) and 



3-dimensional foams (Evans et al. 1995) with re-entrant cell geometries find 
auxetic behaviour. Auxetic materials have many applications, including im- 



pact resistance and sound and vibration damping (Alderson and Alderson 



2007). 



Disordered materials can also be auxetic, but their tensorial description 



is complicated by the inherent inhomogeneities and anisotropics (Horrigan 



et al., 2009). An approximate continuous description can be made using 
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materials exhibit some defects; thus even for highly ordered materials it is 
important to understand how inhomogeneities affect the elastic properties. 

We model a cellular material as a network of elements connected together 
at nodes. For two-dimensional honeycombs, the elements might be beams 
representing the edges of cells. Foams have three dimensional cells; these 
can be open (as in sponges) or closed. However, because of surface tensions 
during manufacture or growth, most of the mass density resides in the edges 



where faces intersect (Gibson et al. 1982). Thus many three dimensional 
foams can still be modeled as a network of beams connected at nodes. We 
assume the network is in equilibrium with no internal forces. We then analyze 
the elastic properties of the network when subject to motions such as hinging 
which make little change to the lengths of the beams. 

The principal mathematical objects of elasticity theory, stress and strain 
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tensors, are locally defined. Most experimental procedures, however, involve 
some sort of non-local averaging. For example, many studies employ the 
engineering strain e^, which for a rectangular material sample gives the ratio 
of the average displacement in the i direction to the length of the region in 
the j direction. 

For a cellular material the microscopic stresses and strains vary strongly 
on the scale of a few cells. Thus an appropriate theoretical description of the 
physics at the cellular scale (or larger) should also involve averaged quantities. 
This paper reports on the properties of a variant of the engineering strain, 
where quantities are averaged over all paths taken by the stresses through 
the material, weighted by the strengths of the paths. This method also 
provides simple interpretations of some basic theorems regarding the mean 
stress tensor. 

Section 2 describes the flow of forces through cellular materials. Sections 3 
and 4 give expressions for mean sizes and displacements in a material sample, 
based on averaging over force paths. Section 5 discusses the calculation of 
the path averaged strain tensor. Example calculations will be presented in 
section 6. In particular, the method of path averaging gives Poisson's ratios 



for regular honeycombs identical to those found in the literature (Masters 



and EvanSj 1996 Gibson and Ashby, 1997 Scarpa et al. 2008). Section 8 



gives conclusions. 

2. Force flows through cellular materials 

The stress tensor determines how forces propagate through a material. 
In the absence of body forces the stress has zero divergence (here repeated 
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Figure 1: The flow of forces in disordered honeycombs. Section |6~3| describes the initial 
setup and boundary conditions for these equilibrium configurations. The left diagrams 
show the directions of F x , while the right diagrams show F y . The size of the arrows is pro- 
portional to the strength of the forces. The equilibrium configurations are initial stressed 
by a small stretch in the vertical direction. Top: A honeycomb with positive Poisson's ratio 
V12 = 0.774, using equation (??). The flow of both F x and F y forces are carried by beams 
under tension, as the directions are rightwards and upwards. Bottom: A honeycomb with 
negative Poisson's ratio z^2 = —0.605. Here the F y forces are carried principally by beams 
under tension, while the F x forces are carried principally by compressed beams. 
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indices are summed and <9j = d/dxi): 

d t cr tJ = 0, J = 1,2,3. (1) 

Letting x be a unit vector in the x direction, we can then construct a diver- 
gence free vector field 

F x = x ■ cr; F xi = a u . (2) 

Similarly, F y = y ■ cr and F z = z • a. The integral curves of F x , F y , and 
F z can be used to track the paths of x, y, and z forces through a medium. 
These curves are consistent with the definition of load path given in |Kelley 



and Elsley (1995). Kelley and Elsley point out that the load path defined 
this way differs from simply following the principal stress direction. This 
difference becomes strong when substantial shear stresses are present. 

Note that F x , F y , and F z can be read off of rows (or columns) of the 
stress tensor. The rows of a tensor are not vectors - but x ■ a is a vector, so 
the latter form needs to be used in coordinate transformations. 

At first sight, one might expect a 180° ambiguity in the direction of a 
load path vector. For example, consider a beam aligned in the x direction 
squeezed from both sides (subject to axial compression). There is no asym- 
metry between left and right; yet the vector F x will point to the left. Here 
the choice of dotting the stress tensor with x rather than — x determines 
the direction of F x . As an is negative for compression, the load flow vector 
F x points in the negative x direction. Under tension F x points to the right. 
Similarly, a vertical beam under tension has F z pointing upwards, etc. 

Consider the forces on a network of interconnecting beams, as in figure 
[TJ Tracking the flow of forces through the beams is analogous to tracking 
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the flow of water through a network of pipes. At each node the x-forces 
(for example) balance; some of the beams connected to a node will have F x 
pointing into the node, others out of the node. The net strength of F x on 
incoming beams equals the net strength of F x on outgoing beams, just as in 
a water network the net flow on inflow pipes equals the net flow on outflow 
pipes. While we can consider the flow of forces in any coordinate system, we 
shall later pay special attention to coordinates aligned along the principal 
axes of the mean stress tensor. 

Suppose the network is enclosed within a boundary surface with outward 
normal n. A subset of the beams cross the boundary. Consider a beam 
labelled a which passes through the boundary at position x(a), with orienta- 
tion along the unit vector r(a). The beam has cross-sectional area A. This 
beam transfers a net x- force $ x (a) through the boundary equal to the flux 
of the load vector F x (a): 



We assume that external forces enter and exit the network only through 
these beams. By force balance, the sum of the external forces is zero. Thus 
the net flux of F x through the boundary vanishes: 



As we have a discrete system, we can write this equation as a sum over the 
boundary beams (see figure [T]). Call the beam a an inflow beam for x forces 
if F x (a) ■ h < 0, or an outflow beam if F x (a) ■ n > 0. For example, if external 




(3) 
(4) 




(5) 
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forces stretch the network in the x direction, so that most of the network is 
in tension, then inflow beams will be on the left, while outflow beams will be 
on the right. Suppose there are N° ut outflow beams, and inflow beams. 
Then by overall force balance 

N m N out 

£$,(a) = £>*(&) = <&,, (6) 

a=l 6=1 

Here the net x- force $ x equals the net inflow (as well as the net outflow). 

The forces entering the system at the inflow beams take multiple paths 
through the network before exiting through the outflow beams. Let $ x (a, b) 
be the net amount of x-force which passes from beam a to one particular 
outflow beam b. Then we have 

N out N in 

$ x (a) = £ $ x (a, by, $ x (6) = $ x(a, b). (7) 

6=1 a=l 

3. The size of a network averaged over force paths 

Consider a particular load path taken by the x-force F x . Various numbers 
can be calculated characterizing the geometry of this path, such as its net 
size in a particular direction. We can then average these numbers over all 



load paths. A well-known theorem (e.g. Landau and Lifshitz (1986), ch. 1.2) 
relates the average stress in a material to a boundary integral over external 
forces. This theorem will be recast in terms of average path sizes and net 
fluxes of the force vectors. 

We will denote an average over F x paths by the symbol ( ) . In this section 
we consider the mean size in the x direction; the next section discusses how 
this mean size changes during a deformation. 
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3.1. The mean path size 

A path from inflow beam a to outflow beam b travels a net distance in 
the x direction 

L x (a,b) = x(b)-x(a). (8) 

We can define the average x distance travelled by the x force by summing 
over all paths the force takes through the network, weighting by the flow 
strength: 

x a=l 6=1 

For a rectangular sample of material aligned with the axes, all paths which 
stretch from the left side to the right side should have the same net length L x 
in the x direction. So why should we seek a path-averaged length like (L x ) x 7 
First, we will need these expressions when relating the mean stress tensor to 
a suitably defined path-averaged strain. Secondly, force paths do not always 
stretch from one side of a rectangular box to the opposite side. For example, 
some x-force paths starting on the left side (carrying a perpendicular external 
load) may end on one of the sides y = constant or z =constant. The x-force 
crossing the boundary at one of these sides connects to an external shear or 
frictional force. Third, quantities like (L x ) x provide a measure of average 
size for irregularly shaped samples. 
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The definition of mean size can be rearranged into a boundary sum: 

= £X>*M)(*(&) -*(<»)) (10) 
0=1 6=1 

N° ut Nl n N° ut 

= J2x(b)J2^> x (a,b)-J2x(a)J2^ x (a,b) (11) 

6=1 a=l a=l 6=1 

N° ut 

= ^x(6)$ x (6)-^x(a)$ x (a). (12) 

6=1 a=l 

We next relate this expression to the mean stress. 

3.2. Boundary integrals 

The mean stress satisfies a useful relation: The stress averaged over an 
area (or volume in 3D) is equal to a boundary integral over external forces 



(Landau and Lifshitz| 1986). Let V be the volume (or area in 2 dimensions) 



o-ij = - I o-ij d 3 x. (13) 



of the sample of material. Then the mean stress tensor is defined as 

1 

V 

As o~ij is divergence free, 

1 f dx~ . 1 



1 
V 



f dx If 

J <7ik dx~ ^ = V J ( (7ikXj ') nk ^ x ^ 
F { ■ h)xj d 2 x. (15) 



Let us look at the 11 component of equation (15). At an inflow rib a 
(integrating over the rib cross-section) 

/(MdWW.) (16) 



while at an outflow rib b 



/(F,.*),d", = *,( W1 ). (17) 
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Thus 

N° ut N' n 



^ = h E ^( fo ) $ x(&) -^E x(a)$ x (a). (18) 



V ^ w w V 

6=1 a=l 



We can write this in terms of the path-averaged (L x ) as in equation (12) 



Similarly 



^=^<^>,; oS=y<^>,. (20) 



4. Path-averaged displacements 

We can treat extensions of a network during a deformation in the same 
way. Suppose that position x changes to x + u. Then the size L x (a, b) of a 
path will change by an amount 6L x (a,b) = Ui(b) — Ui(a). We can average 
this extension over all paths as before: 

N™ N° ut 

( 5L *)* = E E $x ( a > 6 ) 5L *( a > 6 )- ( 21 ) 

x a=l 6=1 



We can rearrange this as in equation (12), then relate it to the boundary 
stress: 



N° N" 



$x(SL x ) x = ^m(6)$ x (6)-^m(a)<l> x (a) (22) 

1 a=l 

(F x • nK d 2 x. (23) 



In general we can define the tensor 
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As the stress is divergence- free, this becomes 

Tij = ^ j> n k a ik Uj d 2 x = ]^j)(Fi- n)uj d 2 x (25) 

= iMMjh- ( 26 ) 

5. The path averaged strain and Poisson's ratio 

We will use our path-averaged extensions and displacements to define a 
mean strain. First, we define the mean strain tensor e in a coordinate system 
aligned with the principal axes of the mean stress tensor. If desired, e can 
be transformed into any other coordinate system. The principal axis system 
provides the most natural context for defining force paths and the fluxes $ x , 
$ y , and $ z employed in the weighted averages. Quantities in the principal 
axis system will be denoted by a *. 

Let x, y, and z be principal axis coordinates. Then we define 

1 = 11 = Wh' ( ] 

and similarly for €2* and €3*. 



Note from equation (19) that (letting a\ = an, etc. 



. g, (5L X ) X * y {SL y ) y $ z (5L Z ) Z 

The general strain tensor, valid in any coordinate system, with these 
quantities as diagonal elements in the principal axis system is 

tij = ((^ife) _1 7fcj + (o^jk^Tki) . (29) 

The Poisson's ratio in the x and y directions constructed from these 
quantities is 

^ = ~- (30) 
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Figure 2: Left: The flow of the x— force vector F x through three connected beams. There 
are inflows along beams 1 and 2, and an outflow along beam 3. Right: The flow of the 
y— force vector F y through the node. There is an inflow along beam 1, and outflows along 
beams 2 and 3. For this particular example, the lengths are proportional to (^1,^2,^3) = 
(27,16,15), and the orientations are (cj>i, <p2, $3) = (— 37r/4, — tt/3, — 7tv/8). The Poisson's 



ratio defined by equation (30) is v\2 = —0.178. The path-averaged horizontal extent is 
$ x (L x ) x = $ x ( pa th a)Aa;(p ath a) + $ x (path p)Ax( path p)- Similarly, in the vertical direction 

*y (Ly) y = ^y(path 7)^J/(path 7) + $y (path 5)Ay( pat h «)• 

6. Examples 

6.1. Three beams connected to one node 

The simplest non-trivial network consists of three beams connected to 



one node (see figure [2]). By force balance, the three beams must lie in a 
plane. Such triply connected nodes form the basic building blocks of many 
two dimensional honeycombs (e.g. hexagonal honeycombs). 

Let the three beams have lengths 2£ 1; 2£ 2 , and 2£ 3 and cross-sections Ai, 
A2, and A3. The half-lengths £1/2, and £3 give the extent of the beams in 
the neighbourhood of the node. The beams have orientations <p2, and 03 
with respect to the x axis. The strength of beam 1 is r 1; etc. Force balance 



13 



requires 



T\A\ COS 0x + T 2 A 2 COS 2 + T3A3 cos 3 = 0; 

TiAxsin^x + r 2 v4 2 sin02 + T3^3sin0 3 = 0. 



(31) 
(32) 



The mean stress tensor is 



Vw = £ x A x o x +£ 2 A 2 o 2 +£ z A 3 o^ 



COS 2 

cos (pi sin < 




6.2. The regular hexagonal honeycomb 

Consider a regular hexagonal honeycomb with beam 1 horizontal (0x — 0) 
and beams 2 and 3 at equal and opposite angles <p 2 = —03 = with respect 
to the x axis. Let £\ = h/2 and £ 2 = £3 = 1/2. To maintain force balance 
at the node, r 2 A 2 = t 3 A 3 = TiAi/{2 cos0). We average the stress over the 
neighbourhood of the node, which includes half of the lengths of each beam. 
Thus 

Tx^lx I h + ^COSi 









(34) 



sm 



COS 1 



The ratio of the diagonal terms is 

oTT _ cos0((/i/T) + cos0) _ sm6((h/£) + sin6) 



V\ 2 



(35) 



cr 2 2 sin cos 2 9 

identical to previous expressions for Poisson's ratio in regular honeycombs 



(Masters and Evans, 1996 Gibson and Ashby, 1997 Scarpa et al. 2008) for 



flexed or hinged ribs. 
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6.3. Disordered honeycombs 

A simple numerical model illustrates the use of path averaging for disor- 
dered two dimensional networks. We start with a regular honeycomb residing 
inside a rectangular region, and perturb the position of each node with a ran- 
dom step taken from a Gaussian distribution. The elements between nodes 
are assumed to carry a linear (Hooke's law) force. The outermost nodes 
are then moved outwards or inwards slightly, and the equations of motion 
are solved with heavy damping until the system settles into an equilibrium 
state. The mean stress tensor can now be calculated, giving a prediction of 
the Poisson's ratio of any small motions to or away from this equilibrium. 
To obtain the displacement, we move the top and bottom boundary nodes a 
tiny amount, while letting the nodes on the left and right boundaries move 
freely. 

To illustrate for a reentrant honeycomb, we began with the cells in a 
regular hexagon pattern with i = h/2, <p = 0.657T, and U21 = —0.884, and 
perturbed with random steps (rms O.l/i) (see bottom of figure [T]). A stretch 
in the vertical direction of O.lh was followed by a relaxation to an equilib- 
rium. The equilibrium configuration had V21 = —0.605. The inclination of 
the principal axis frame was negligible: only 0.02 degrees. The individual 



nodes can be analyzed individually, as in section 6.1 Here the average v of 
the individual nodes was —0.566 in their principal axis systems, not much 
different from the total configuration. However, the principal axes of individ- 
ual nodes were often misaligned with respect to the total: in the xy frame, 
Poisson's ratios varied from -4.28 to -0.45. The Poisson's ratio calculated 
without weighting by path strength was 0.595, only two percent off from the 
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path-averaged ratio. 

For the top example, the cells started as regular hexagons stretched in 
the horizontal direction, with h = 1.5£, <ft — tt/3, and v 2 \ = 0.75, before 
the random perturbation. The Poisson's ratio became v 2 i = 0.774 after a 
random perturbation. The inclination of the principal axis frame was 15 
degrees. The average z/ 2 i of the individual nodes was 0.685. The Poisson's 
ratio calculated without weighting by path strength was 0.903. 



7. Conclusions 

The elastic properties of a disordered cellular material depend on the 
size and geometry of the sample being measured. To describe these elastic 
properties, we need to employ some method of averaging in order to obtain 
effective elastic moduli. This paper has examines the possibility of averag- 
ing over all paths taken by the forces through the material. This method 
emphasizes the paths which carry the most stress, and hence have the most 
influence on the mechanics of the material. Also, in many cases it may give 
effective moduli close to those calculated by cruder averaging. 

This result has immediate implications for engineers wishing to design an 
auxetic material. Consider a honeycomb stretched in the y direction, so that 
the mean 022 < 0. Then for v\ 2 to be negative, we must have a positive 0T1, as 
in figure [Tj Thus the solution to the force balance equations must provide a 
mixture of elements where some are compressed and some are under tension. 
Here the y forces must be primarily carried by the elements under tension, 
and the x forces carried by the compressed elements. 



A reentrant honeycomb geometry (Masters and Evans, 1996) accom 
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plishes this naturally: while the y force F y must ultimately flow from bottom 
to top, the paths zig-zag both upwards and downwards on the way The 
downwards sections of the paths are under compression; these contribute to 
the negative flow of x forces F x . Not all auxetic networks have an obvi- 



ous reentrant geometry (Horrigan et al. , 2009); an analysis based on load 



paths and the mean stress tensor may help to understand auxetic behavior 
of disordered networks in general. 

Furthermore, there may be a hierarchy of structures contributing to aux- 



etic behavior (Taylor et al. , 2011). The methods proposed here can be mea- 
sured for any subvolume of a material, rectangular or irregular in shape. This 
flexibility allows analysis of elastic properties on many length scales. 
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